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Abstract. Auger results on clustering of £ 60 EeV ultra- high energy cosmic ray (UHECR) 
ions and the interpretation of the 7-ray spectra of TeV blazars are connected by effects from 
the extragalactic background light (EBL). The EBL acts as an obscuring medium for 7 rays 
and a reprocessing medium for UHECR ions and protons, causing the GZK cutoff. The study 
of the physics underlying the coincidence between the GZK energy and the clustering energy 
of UHECR ions favors a composition of £ 60 EeV UHECRs in CNO group nucleons. This has 
interesting implications for the sources of UHECRs. We also comment on the Auger analysis. 



1. Introduction 

The recent reports [U [2] from the Pierre Auger Observatory collaboration of high-significance 
clustering of UHECRs along the supergalactic plane has opened multi-messenger astronomy. 
This advance is closely connected to progress in 7-ray and neutrino astronomy. Here we focus 
on new techniques to unravel the problem of UHECR origin in light of the Auger results. This 
topic is developed in more detail in my Merida ICRC contribution [3] and in my book with 
Govind Menon [JJ. 

2. Charged Particle Astronomy 

We begin by reviewing the recent major results from Auger, including the discovery of the 
clustering of the arrival directions of UHECRs along the supergalactic plane [I] . This discovery 
was anticipated in the study by Stanev et al. [5j, and we return to this paper in the discussion. 

2.1. Pierre Auger Observatory 

Auger is ~ 3000 km 2 in area, about the size of Rhode Island. It consists of 1600 shower detectors 
spaced 1.5 km apart, as well as four telescope enclosures each housing 6 telescopes to map Ni 
fluorescence. The directional point spread function of an UHECR is typically ^ 1°. The original 
analysis published in Science pQ treated a data set from 1 Jan 2004 to 31 August 2007 consisting 
of 81 events with total energy E > 40 EeV. 



2.2. Clustering Result 

A probability statistic P corrected for exposure is constructed from the nearest-neighbor angular 
separation ip between the arrival direction of an UHECR with energy E and the direction to 
AGN in the Veron-Cetty & Veron (VCV) catalog [6j. This catalog has 694 active galaxies with 
z < 0.024 or distance d ^ 100 Mpc (the total catalog has over 10 5 entries) but is incomplete, 



especially around the plane of our Galaxy, and is also biased by different sample selection 
techniques. 

In the Auger analysis [lj, P was minimized for ip = 3.1°, threshold energy Ethr = 56 EeV, 
and z max = 0.018 (d max = 75 Mpc), containing 27 events (the two highest energy events were 
90 and 148 EeV). Twelve events correlate within 3.1° of the selected d < 75 Mpc AGNs, and 
another three within the vicinity of one of these nearby AGN. By comparison, 3.2 events would 
be expected for an isotropic UHECR flux, which it clearly is not. Lack of correlation of £ 60 
EeV UHECRs with the Galactic plane rules out Galactic (pulsar, stellar mass black hole) and 
Galactic halo (including halo dark matter annihilation) models, leaving extragalactic models 
viable, at least those that follow the matter distribution traced by the supergalactic plane like 
the nearby AGN in the VCV catalog. 

The recent detailed Auger analysis [2] confirms these results, furthermore pointing out that 
the distance cutoff could range from ~ 50 - 100 Mpc and that ip £ 6°. If one accepts the validity 
of these results, then UHECRs with E ^ 60 EeV originate from AGNs within 75(±25) Mpc and 
are deflected in their travels by at most ~ 3° - 6° by either the Galactic magnetic field [2] or 
the intergalactic magnetic field [3]. 

The lack of a quasi-isotropic background UHECR flux above E > E c \ = 60 EeV is most 
simply explained if more distant UHECRs fail to arrive here at Earth. Because we know that 
there are sources of UHECRs with energy up to at least 300 EeV [7], it follows that unless 
we were in a unusual over-density of UHECR sources within ~ 100 Mpc and an under-density 
within ~ 100 - 300 Mpc or more, then energy losses of UHECRs accelerated by and injected 
from more sources more distant than ~ 50 - 100 Mpc are the cause of the Auger anisotropy. 

The clustering energy E c \ = 60 EeV stands out in this analysis, separating E ^ E c [ UHECRs 
formed mainly in the large scale quasi-homogeneous, isotropic universe at distance scales d £ 50 
- 100 Mpc from the UHECRs formed in the clumpy structured universe at local, d ^ 100 Mpc 
scales. 

2.3. GZK Cutoff 

A high-significance steepening in the UHECR spectrum at the Greisen-Zatsepin-Kuzmin (GZK) 
[SI E] energy E GZK ^ 10 19 - 6 eV ^ 4 x 10 19 eV was reported at the 2007 Merida ICRC based on 
observations taken with the Auger Observatory [10] (see Figure [1]) , and earlier in 2007 by the 
HiRes collaboration p] at E GZK 10 19 - 8 eV = 6 x 10 19 eVW This observation seems to confirm 
the prediction [8j [9] that interactions of UHECRs with CMBR photons will cause a break in the 
UHECR spectral intensity near 10 20 eV. The coincidence E GZK ^ E cl seems likely to originate 
from underlying physics. 

24. UHECR Composition 

The composition of UHECRs measured from ~ 4.5 x 10 17 eV up to 4 x 10 19 eV with Auger is 
neither dominant p nor Fe composition [12], and is not trending toward a lighter composition 
at the upper energy range, contrary to pre- Auger indications |13| . If the results of Ref. [12] can 
be extrapolated to 60 EeV, then ions must be included in the analysis unless the UHECR 
composition abruptly shifts to proton-dominated at E ^ E c i, as would be the case in neutral 
beam models of UHECRs [H]. 

3. GZK Energy 

The energy of the GZK cutoff is not precisely defined. Returning to the classic paper by Greisen 
[8J for inspiration, let us work through the estimates in light of our new knowledge. 

1 The HiRes team also reports a significant dip structure at Edip — 10 18 ' 6 eV. 
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Figure 1. UHECR 
spectrum, plotted in the 
form of E^I, where I 
is the number intensity, 
for AGASA (filled di- 
amonds), HiRes (open 
symbols), and Auger 
(filled circles). 



3.1. GZK Energy: Analytic Analysis 

UHECR protons with energy E p = jpiripC 2 = 10 20 £ ; 20 eV will undergo strong photopion 
losses when the threshold condition j p €cmbr ~ m^/rrie = 400 is satisfied, where ecmbr — 
2.70 x 2.72 K/5.9 x 10 9 K = 10 -9 is the dimensionless mean photon energy, in units of the 
electron rest mass, of the CMBR. Strong photopion losses due to interactions with the peak 
energy photons of the CMB occur when 7 p }t 4 x 10 11 . The mean free path (MFP) for energy 
loss of a proton with energy E20 ~ 4 is A^- = {l/ncMBR^<j>w) — 1/(400 cm -3 • 70/zb) = 10 Mpc. 
From a derivation using a step- function approximation for the product of photopion cross section 
and inelasticity [13], the mean- free-path for photopion energy losses of UHECR protons with 
CMBR photons is given by [3] 

^ 13.7exp[4/E 20 ] 
X ^ {E2o) = [1 + 4/E 20 ] MPC • (1) 

The values of A^E^o) for E20 = 0.6, 0.8, 1, and 2 are 1400, 340, 150, and 35 Mpc, respectively. If 
the UHECR with energies between 60 and 100 EeV were protons, then they could originate over 
a large distance scale extending to 1 Gpc. This simple estimate suggests that the the UHECRs 
cannot be protons because the Auger events originate from within 100 Mpc whereas the energy 
loss MFPs of UHECR protons at these energies are much greater. 

Now we estimate the GZK MFPs of ions against photodisintegration on CMBR photons. 
Note that by contrast with UHECR protons, photopion losses of UHECR nucleons are not 
effective until r y(e)cMB ~ m n /m e or £20 ~ + z), where A is the nucleon mass number. 

The pion-forming A resonance peak energy is = m e c 2 e r ~ 300 MeV, a much higher energy 
than the giant dipole resonance (GDR) at ~ 20 MeV, so photodisintegration breakdown usually 
dominates photopion losses for ions, especially for large A. The effective photodisintegration 
energy-loss rate for ions on the CMBR results from both GDR excitation with the emission of 
one and two nucleons, and from multi-nucleon breakup processes. Because the reaction is a 
threshold process and UHECRs with E20 <C 3A/(1 + z) interact only with the Wien portion 
of the CMBR radiation field, the GDR is more important than multi-nucleon excitations for 
£20 ~ 1 UHECRs. 

In this asymptote, we use the 5-function approximation for the photodisintegration cross 
section of nucleus A, namely 

o"A(er) — ^o.iAcflfl^er - e r> A) , (2) 



PS HUE], where 



a ,A = 1.454 mb , A GDR ^ 15.6 , e r , A = 83.5A -0 ' 21 , 

and e r = — /i) is the invariant dimensionless photon energy in the ion's rest frame. In the 
5- function approximation for the GDR, eq. ([2]), the derived inverse of the effective energy loss 
time scale of an UHECR ion due to photodisintegration processes with photons of the CMBR is 

t E (E, A) = -— ^ In ([1 - exp{-w r} A)] ) 

,= k A A 1 - 79 (l + z) , _ 1 

3.2 x 10~ 15 — „ V ; In , s , 3 

E 2 L l-exp (-w r ,A 



where 



e r ,A 



0.834° 



79 



2 7 G £ 20 (l + z)" 

The effective path length for photodisintegration in the <5-function approximation for the GDR 
is therefore 

x WP A s . (w as 3.Qgf exp(0.834 - 79 /[(l + z)E 20 ]) 

A B = A S (.E,A) = ct E (E,A) -> fc A A 1 - 79 (l + ^) ^ 

in the limit w r> A » 1, that is, ^20 < 5(yl/10) - 79 . 

When a single proton or neutron is ejected, then a fraction A^ 1 of the original energy E 
is lost to the original nucleon, and for the ejection of two nucleons, a fraction 2 J A of energy 
is lost. For multi-nucleon injection, an average fractional energy loss k A /A is used [T7j, where 
fc,4 = 1.2/A,3.6/A, and 4.349/A for A = 4, 10 < ^4 < 22, and 23 < A < 56, respectively. The 
photodisintegration energy-loss MFPs have only a generalized meaning in terms of the mean 
distance for an UHECR nucleon to be broken up into mostly lower energy protons and neutrons 
and a nucleon with A about half of the original nucleonic mass. Photodisintegration of Fe, for 
example, leads to significant fraction of A £ 15 nucleon secondaries [T7). 

Equating eq. (H|) with an energy-loss distance A = IOOA100 for E = 60E , 6o EeV gives 
exp(1.38A a79 /£ 60 )/[/cAA L79 ] = 333Aioo^Io- li is easil y verified that all nuclei with A > 14 
have energy-loss mean free paths at Eqq ~ 1 much longer than 100 Mpc. On the basis of 
this analytic treatment, one might conclude UHECRs could originate only from nucleons with 
4 < A < 14. 



3.2. GZK Energy: Numerical Analysis 

These results are not conclusive unless additional radiation fields, including IR and optical 
features from stellar and reprocessed radiation, and energy-loss processes of photopair production 
and universal expansion are taken into account. Figure [2] shows the energy density of the EBL 
and cosmic rays, and the predicted upper limit to the cosmogenic neutrino energy density |21j . 



3.2.1. Extragalactic Background Light The mean intensity of light in intergalactic space is 
referred to as the EBL. The IR and optical EBL is decomposed at low redshift {z £ 0.25) into a 
dust component and a stellar component as shown in Figure [3l These fits are assumed to span 
the likely range of the diffuse radiation fields between ~ 1 and a few hundred /im at low z. The 
modified blackbody spectral energy density is written in the form 

eu(e) = u — = m e c 2 e 2 n ph (e) , (5) 

exp(w) — 1 
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Figure 2. Energy den- 
sity in intergalactic space 
of various radiations, in- 
cluding the CMBR, the 
IR and optical, X-ray, 7 
ray, extragalactic cosmic 
ray, and the predicted 
maximum energy density 
of cosmogenic neutrinos. 
Also shown is the near- 
Earth energy density of 
cosmic rays; the transi- 
tion energy between the 
galactic and extragalac- 
tic component is uncer- 
tain. 



Table 1. Properties of the Dust and Two Stellar Components 



Component 


T(K) 


(10 -14 ergs cm -3 ) 


k 


CMBR 


2.72 


6.38 


4 


Dust (Hi) 


62 


0.819 


3.8 


Dust (Low) 


31 


0.273 


3.8 


Star 1 (Hi) 


7100 


1.1 


2.0 


Star 1 (Lo) 


7100 


0.55 


2.0 


Star 2 


16,600 


0.5 


3.0 



where w = e/Q. The stellar component is decomposed into the sum of two modified blackbodies, 
with the higher-temperature stellar component fixed. The HI EBL has the more intense dust 
and stellar fields, and the LO EBL has the weaker dust and stellar fields. Two other EBLs are 
considered here: a high dust/low stellar field and a low dust/high stellar field. Parameters of 
the CMBR and modified blackbodies are given in Table 1. 

The intensity of the diffuse and unresolved 7-ray EBL from EGRET observations and analysis 
|23[ [2"4"j is shown in Figure [H It is decomposed into contributions from blazars, separated into 
FSRQs and BL Lacs, and the superposed intensity from numerous faint sources where the 7 rays 
are formed primarily by interactions of cosmic rays accelerated by supernovae in star-forming 
galaxies and by structure formation shocks in clusters of galaxies (reviewed in [26\ |2"7]). The 
particle and radiation fields are connected because GZK effects of UHECRs with the EBL form 
cascade 7 rays that contribute a truly diffuse component to the 7-ray EBL [25]. The intensity 
of this component is sensitive to the transition energy from galactic to extragalactic UHECRs 
and the formation rate of UHECR sources through time. 

The analytic result for the energy-loss MFP of Fe with A = 56 on CMB photons is given 
by eq. dH) and labeled "analytic" in Figure [5j The accuracy is tolerable, but a more detailed 
calculation, also shown in Figure [5l is needed to draw more reliable conclusions. The numerical 
calculations for Fe treat photopion and photopair losses with the EBL as well as from universal 
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Figure 3. Ranges of EBL used 
in analysis. Data from Ref. |22j . 
The HI (LO) EBL refers to a high- 
(low-) intensity dust component 
and a high- (low-) intensity stellar 
component. 
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Figure 5. Effective energy-loss 
mean-free-path of Fe in LO and HI 
EBLs, with the effective photodis- 
integration energy-loss MFP calcu- 
lated according to the prescription 
of Ref. p3]. 
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expansion (the Hubble radius is ~ 4160 Mpc at the present epoch), with the photodisintegration, 
photopair, and photopion components shown separately for the HI EBL. The total energy loss 
MFP is also shown for the LO EBL. (See also [281 129j.) The distances from which we should 
expect to detect 60 - 100 EeV Fe range from « 200 Mpc to 150 Mpc. 

Fig. [6] shows the energy-loss MFP of UHECR protons in the LO and HI EBLs. The MFPs, 
which include photopion, photopair and expansion losses, range for the HI EBL from 200 to 100 
Mpc for energies ranging from 60 to 100 EeV, respectively (similar to the values for UHECR 
Fe). Unless there was a deficit of UHECR sources between ~ 100 - 200 Mpc, or an enhanced 
pathlength due to strong, ^ 10~ 10 G intergalactic magnetic fields, it is difficult to account for 
unless UHECR p and Fe comprise only a small fraction of the ^ 60 x 10 19 eV UHECRs detected 
with Auger. 

How else to resolve the quandary that the UHECR proton and Fe MFPs are too long to be 
commensurate with the Auger data, so that the UHECRs cannot be protons or Fe? Possibly the 
Auger detector is not properly calibrated, though it is unlikely to be off by more than ~ 20% 
compared to the factor of « 2 needed. For instance, Berezinsky [30], using the dip energy seen 
most clearly in an E 3 I plot (see Fig. [T| as a fiducial, and multiplying by an energy calibration 
factor 1.2, 0.75, 0.83, and 0.625 for Auger, AGASA, Akeno and Yakutsk data, respectively, puts 
the various data sets in agreement with each other. Watson [33] argues that there are ~ 19% and 
11% errors in the Auger energy calibration at 10 17 and 10 20 eV, respectively, so the correction 
factor for Auger is 1.2. If this is true, then only a few percent of the ^ 6 x 10 19 eV UHECRs 
could be protons. 

Let us provisionally accept that the Auger clustering results and the calculations presented 
here require that ^ 60 EeV UHECRs are depleted in p and Fe (a full Monte Carlo simulation 
is required to remove remaining uncertainties). Some models, in particular, the neutral beam 
model that Armen Atoyan and I developed (for blazars, see [H], and for GRBs, see [31]) were 
based on the acceleration of protons to ^ 10 20 /r eV energies in the shocked fluid frame, where 
r is the flow Lorentz factor. Photopion processes make an escaping neutron and 7-ray beam. 
If protons do not make up a large fraction of the ^ 10 19 ' 6 eV UHECRs, then this model no 
longer has to contend with the difficulty of accelerating protons to such high energies. On the 
other hand, this or any other model must be able to accelerate ions to ultra-high energies and 
get them from the accelerator into intergalactic space with the requisite emissivity. 

Which ions? Again accepting the approximate integral meaning of an energy-loss MFP, 
Figures [7| and M show that UHECR ions up to Ne (A = 20) or at most Mg (A = 24) would 
effectively lose their energy within ~ 100 Mpc. An abundance of UHECRs in light ions gives the 
most economical explanation of the Auger results. If correct, then composition analyses of 60 - 
100 EeV UHECR data will give an average atomic mass number (A) ^ 24 ((In A) ^ 3.2), and 
probably in the range (A) ^ 20. This prediction can be compared with Fig. 10 in Watson |33| . 
giving the measured energy dependence of (lnvl) using the the QGSJETII-03 model. According 
to our calculations, the anisotropy of = 60 EeV UHECRs with d max ^ 50 - 100 Mpc could only 
originate from muclei with 2 ^, (ln^l) ;$ 3 (the lower limit is not firm, and depends on source 
type). 

Let us return to the question of the EBL. The dust EBL primarily attenuates ^ 4 TeV 
photons and the stellar EBL primarily attenuates ^, 4 TeV. The LO and HI EBLs can be used 
to calculate r 77 up to redshifts z « 0.25 where cosmic evolutionary processes start to change the 
EBL energy density compared to the present. The exponential attenuation factor exp(r 77 ) as 
a function of photon energy E for different redshifts is plotted in Figure [9l We can make a cut 
defined by exp(r 77 ) = 10 (or any other number, for that matter) to define a relationship between 
redshift and photon energy (the relation defined by r 77 (z,E) = 1 is sometimes referred to as 
the "Fazio-Stecker relation [M]"; see also [35]) Here we make the cut defined by exp(r 77 ) = 10 
to give the obscuration energy E osc as a function of obscuration redshift z osc , shown in Fig. 
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Figure 7. As in Fig. ([5]), for 
UHECR N. 



Figure 8. As in Fig. (jSJ), for 
UHECR Ne. 



10* 



z = 0.025 to 
z = 0.25 

in steps of 0.025 



HI EBL 




Figure 9. The factor exp(r 77 ) at 
different redshifts for the LO and 
HI EBLs shown in Fig. | as a 
function of photon energy E. 
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Figure 10. Relation between 
the obscuration energy E osc and 
obscuration redshift z osc defined by 
exp(r 77 ) = 10 for various high and 
low dust and stellar EBLs. 



[TU1 This figure is not however very useful for comparison with observational data, because the 
highest measured energies of photons depend upon telescope sensitivity, and variable source 
spectral brightness. Nevertheless, the dust component forms an opaque screen to photons with 
E > 7 TeV and > 15 TeV for the high and low dust EBL at the redshift z « 0.03 of Mrk 421 
and Mrk 501. Significant detection of w 8 TeV from Mrk 421 [361 EI] and Mrk 501 [38], and of 
> 10 TeV photons from 1ES 0229+200 [39] (z = 0.14) suggests that the EBL intensity at ~ 10 
- 100 fim is smaller than given by the high dust EBL. 



We now ask why the sources of UHECRs should be enriched in medium mass (C - Ne) nuclei. 
Refs. [15} 116] discusses various scenarios for accelerating UHECRs with significant ion content in 
GRB environments, including capturing supernova shell material by a jetted relativistic wind, 
external shock capture of circumburst material, and hypernovae [40J. We favor the hypothesis 
|41j that an external shock is important during the prompt and early afterglow phases of a GRB. 
GRBs formed by core collapse of C or O Wolf-Rayet stellar progenitors would be surrounded 
by a wind enriched in these light elements that could be captured and accelerated to ultra-high 
energies. This could be accomplished by 1 st - or 2 nd -order Fermi processes at the shock formed 
by the interaction of a relativistic wind with the surrounding medium. 

For blazars and radio galaxies that propel intermittent winds to form internal shocks, the 
composition of UHECRs would depend on the makeup of the matter accreting onto supermassive 
black holes. If near-Solar metallicities occur [32], then proton acceleration would have to be 
suppressed, perhaps by a gyroresonant mechanism, and the medium mass ions would have 
to be processed through the black hole environment without breakup due to spallation or 
photodisintegration. 

4. Discussion 

The clustering result reported |T] by the Auger Collaboration is a major discovery, anticipated by 
Stanev, Biermann, Lloyd-Evans, Rachen, and Watson [5] in 1995. In the earlier study, Northern 
hemispheric data from Haverah Park, AGASA, Volcano Ranch, and Yakutsk with 73 (27), 25 
(7), 13 (5), and 32 (3) events above 20 (40) EeV, respectively, and with zenith angles smaller 
than 45° were analyzed. The average and rms angular distances of the arrival directions of 
UHECRs from the supergalactic plane were calculated, after correcting for exposure. Compared 
to an isotropic source flux, a statistical enhancement towards the supergalactic plane at the 2.5 
- 2.8 a level was found for arrival directions of events with E > 40 EeV. The supergalactic 
plane runs through the Virgo Cluster at ~ 20 Mpc, and contains an assortment of radio glaxies 
such as M87, Cen A and NGC 315, and the starburst galaxies M82 and NGC 253. Stanev et 
al. [3] argue that their analysis supports the hypothesis that radio galaxies are the sources of 
UHECRs accelerted, for exmple, in the inner jets or at the knots or hot spots of radio galaxies. 
Supporting this contention is the analysis of Shaver and Pierre [12] showing harder integral 
source count distributions of extragalactic radio sources at large spectral fluxes toward the 
supergalactic plane in the 408 MHz 1 Jy southern hemisphere Molonglo catalog, and clustering 
of sources in the 2.7 GHz 2 Jy Wall and Peacock catalog [H] towards the supergalactic plane. 
This clustering conforms with the association of UHECR arrival directions with radio sources 
in the supergalactic plane. 

By comparison, the Auger team [H12] uses the VCV catalog [6j consisting of 85,221 quasars, 
1122 BL Lac objects and 21,737 active galaxies (including 9628 Seyfert Is). Quoting Veron- 
Cetty and Veron, "This [VCV] catalogue should not be used for any statistical analysis as it 
is not complete in any sense, except that it is, we hope, a complete survey of the literature." 
The associations between the arrival directions of UHECRs and directions to nearby AGN in 
the VCV catalog may only reflect that both the AGNs in the VCV catalog and the sources 
of UHECRs within ~ 50 - 100 Mpc are preferentially found toward the supergalactic plane, 
and that the large number of sources in the VCV catalog gives an apparent stronger underlying 
association than would be found if a smaller catalog, such as the Wall and Peacock catalog, were 
used. Associations made with a mixed catalog are more likely to be spurious, and the removal 
of statistical biases in such analyses are challenging. Correlating UHECR arrival directions with 
flux-limited surveys of sources with relatively complete redshift information is the only reliable 
way of identifying the source population(s) of UHECRs in the absence of high-energy neutrino 
detection from these sources or discovery of anomalous 7-ray signatures in radio galaxies and 
GRBs with GLAST. 



A concern expressed recently [45} 06] is that although some of the UHECRs observed with 
Auger are associated with the direction toward Cen A, none seem to originate from the Virgo 
cluster center of our supergalaxy in the direction of M87. The lower exposure of Auger toward 
M87, by a factor of 3 compared the exposure to Cen A [2], and the Northern Hemisphere 
detection of UHECRs that could be associated with Virgo [5j 01] could ameliorate this issue. 
Ref. [46] suggests that if all UHECRs are C, then the short energy-loss MFP of C would mask 
M87 and explain preferential detection of Cen A. But our calculations show that the effective 
energy loss MFP of UHECR C with E = 60 EeV is ps 50 Mpc, within which M87, at a distance 
of ps 18 Mpc, falls. The MFP of UHECR He at 60 - 100 EeV is ~ 3 - 5 Mpc, and so could 
restrict arrival directions to Cen A. 

The dip energy, viewed as a consequence of photopair effects on UHECR protons [30J either 
from GRBs [37] or AGNs [3H], could involve photodisintegration and photopair energy losses 
of ions interacting with photons of the EBL. There is also the possibility, again, that there is 
a deficit of UHECR sources between ps 100 and ps 300 Mpc. Whether this is true depends on 
source type, in particular whether UHECRs originate from specific classes of GRBs or AGNs. 

In summary, the Auger team has opened charged particle astronomy by correlating the 
anisotropic arrival directions of ^ 60 EeV UHECRs with AGNs located within ps 75 Mpc. 
Calculations presented here of the effective energydoss MFP support an origin of UHECR ions 
in light nuclei with A ;$ 24. The conclusion following from this analysis is that detection of 
UHECR clustering is most easily understood if most UHECRs are mainly light nucleons. The 
identification of the actual underlying source population must use complete (at least at high, 
|6| ^ 10°, galactic latitude) fluxdimited catalogs. 
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